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In this paper we calculate additional contributions to that part of the non-Gaussianity of the 
primordial curvature perturbation which come from the three-point correlator of the field pertur- 
bations. We estimate this contribution in the following models for the origin of ^: single-component 
inflation, multi-component chaotic inflation, a two-component "hybrid" inflationary model, and 
the curvaton scenario. In all of these models, the additional contributions to the primordial non- 
gaussianity considered here are too small to ever be detected. 
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I. INTRODUCTION 



There is substantial interest in an observable deviation from a Gaussian distribution of the primordial curvature 
perturbation C [1 II ijl II IE II S i, E 111 III 111 Q 111 

From observation, the nongaussian contribution of C is 
expected to be small [ig- Nevertheless the current bound still includes a wide range of models for the origin of (, and 
therefore a future nongaussian detection or a tighter bound on the nongaussianity of C, could rule out a considerable 
number of models. 

The level of nongaussianiaty in the curvature perturbation can be parameterised by /nl, which specifies the three- 
point correlator of (. Here we apply the 6N formalism 0, 0| to calculate new contributions to /nl- While the 
relevant scales are outside the horizon, the curvature perturbation according to the SN formalism, is given by H 

C(x,t) 6N{p{t),cl>,{^)) (1) 

= N^{t)H^i^) + ^ E N,3{t)SM^)S^,{^) + ■■■ , (2) 

i ij 

where iV(p(t), 0i(x)) is the number of e-foldings from an initial flat slice, on which the field values are evaluated, 
to an uniform energy density slice at a time t after has settled down to the time-independent value. Ni and Nij 
are the partial derivatives, Ni = dN/d(j)^ and Nij — d'^N/d(j)^d4>^ , resulting from the Taylor expansion on the field 
perturbations. From now on, we consider only the first two terms of the expansion in Eq.(|2l 
In the case of Gaussian field perturbations /nl is given by 



Seery and Lidsey [g have calculated the three-point correlator of the field perturbations. To leading order in the 
field perturbations, the three-point correlator is 

(Ck,Ck.Ck3> C E N.N.NkiS^blMMs) ' (4) 

ijk 



which gives the following additional term 



J:^,kN^NJNkf"' (fcl,fe,fc3) 



A/jvL = ''^ \: . (5) 

[j:,N^]y^vi^' 

The /*J'^'s are dimensionless functions proportional to the three-point correlator of the field perturbations. This term 
adds a contribution to /nl too small to ever be observable, |A/nl| < 0.044 Jdj. 

In this paper we calculate the next order contribution to /nl- Contrary to the Afj^^ case, we are unable to give 
an estimate for an arbitrary potential. Therefore we calculate its contribution to the nongaussianity of ^ in several 
specific cases: single field infiation, multicomponent chaotic inflation, the two-component inflationary model of Ref. 
[l9l. and in the curvaton scenario. 



II. HIGHER ORDER CORRECTIONS TO /nl 



The observed distribution of the curvature perturbation ^ in the CMB is almost Gaussian. The power spectrum in 
the SN formalism is then approximately given by 

^ ^ i 

The three-point correlator of or its bispectrum defined by (CkiCk2Ck3) = (27r)'^i?f S^Cki + k2 + ka), is the lowest 
order signature of non-Gaussianity. To measure a possible small deviation from gaussianity in ( we introduce the Jnl 
function, defined as Q 

(fci, fe, ^3) = ^/nl [Pc ih) Pc (fca) + cyclic] , (7) 

where P^{k) = 2tt'^'P c, / . 

We consider the contribution of the three-point correlator of the field perturbations to the following term^ 



(Ck,Ck.Ck3> C i ^ N,N,,Nuk{54>l, {6<j,^)l^ {5^%^) + (k, ^ k2) + (fci ^ fca) , (8) 

ijk 

where ki ^ k2 and ki ^ ks represent the corresponding permutations, and denotes the convolution product of 

the field perturbations: 

' d^qS(f>qS<l)w-q. (9) 



a2\ 



The correlator in the field perturbations has the following generic form 

(SKSKML-^Ms^^L-^^s) = ('50L/0^q.>(^0L-q.<'5C-q3> + Comb, (10) 

where Comb refers to all possible combinations of the 5 elements in groups of two, which in total number 10. In the 
equation above, the two point function is defined by 



fki' 



(27r)'(5^-'"P5.. (fci)<5(3)(ki-fk2) , (11) 



where we consider a scale invariant spectrum 'Ps^i{k) = {k^ /2Tr'^)Ps^i = {HI 2tt)'^ . On the other hand, the three-point 
function obtained from second order cosmological perturbation theory is 



with AI123 given by 

7Wi23 = M (fci,fc2,A:3) = ^ 



-3^-^(fci+2fc3) + ^fc?-fcifc2' 



(13) 



where kt denotes the sum of the modulus of the three momenta kt — ki + k2 + ks. In Eq. p2|) the sum a' is over all 
the permutations of the indices i, j and k, and at the same time over their respective momenta fci, ^2 and k^. Mp is 
the reduced Plank mass: Mp = 2.436 x 10^^ GeV. 

Of the 10 combinations in Eq. pO|l only 8 are non-vasnishing, since neither k2 nor k3 can be zero. The 8 terms can 
be arranged into two groups of 4 terms, each group containing terms that give the same contribution. These can be 
represented schematically as two different diagrams, shown in Fig.l^. The other two contributions in Eq. ((SJ arising 



^ Note that we have diagonahzed the Nij in the field perturbations basis to simplify the calculation. The connected part of the five point 
function is not written explicitely. 



FIG. 1: Here we show the two possible diagrams. There are three external lines corresponding to the three momenta fci, ^2 
and fcs. The parallel lines represent the two-point correlator, whereas the crossed circle indicates the three-point correlator. 
The diagrams have been drawn using JaxoDraw [2nl |. 

from permuting the ki momentum with ^2 and k^, result in the same two diagrams. Similar diagrams are given in 
Ref. j21|]. 

In calculating the contribution to the three-point function of ( resulting from the first diagram, we find the following 
integral over the internal momentum: 



^123 — 



i3 



■M(fci,|k2-KqUk3-q|) 
'3 fc3|k2 + q|3|k3-q|3 



(14) 



and similar integrals corresponding to the permutations in the argument of the M123 function. The momenta de- 
pendance of the integrand can be read off from the location of the three-point correlator in Fig.(^|. In the equation 
above L^^ indicates that the integrand is set equal to zero in a sphere of radius around each singularity. The 
integral then is approximately given at each singularity by 



'123 



47r InKL 



X(fci,|k2-f q|,|k3-q|) 



q|3|k3-q|3 /q=o,-k.ork3 



(15) 



where the function in brackets is evaluated at the singularity, and K = min{/cj}. The logarithm in the expression 
above can be taken to be of order 1 . It can be shown that each permutation of the M. 123 function only diverges at 
q = 0, and only at one of the other two singularities if at any, since the integrand becomes in at least one of them. 
Therefore we are able to give a momenta dependent result for each integral. 

The total contribution to the the Jnl function from the first diagram, using the identity Ni(p'^ — —H as in [TEj l. 
gives 



M2 



The momenta dependent term has a maximum absolute value for k = ki = k2 = k^ given by [TH 



E.M{k,k,k) 



11 



(16) 



(17) 



where the sum over a denotes the sum over the permutations of the three momenta. In the following sections, we 
calculate the value of A/^]^ for the maximum of the momenta dependent term. 

The integrals for the second diagram are evaluated in a similar way to the integrals of the first diagram, but this 
time there are more cancellations at the singularities. The total contribution to f^^ function from the second diagram 
is 



- ^2 1^2^ 



(18) 



In the following sections, we evaluate these higher order contributions in Eqs. (|16|l and H18() for several specific 
cases. 



III. SINGLE-COMPONENT INFLATION 



In single field inflation the number of e-foldings N{t,to) is given by 



1 /"Po V 



(19) 



In the 5N formalism we evaluate the partial derivatives of N respect of the field values on the initial flat slice. Then 

1 



2eM„ 



where e = Mp/2 {V^/VY and 77 = Mp(y^^/V) are the slow-roU parameters. 
Defining A/^^^'* = A/^]^ + A/^]^, we have the following contribution to / 

t 

4 

This is much less than the leading order contribution calculated in 



NL- 



(20) 
(21) 



(22) 



IV. MULTI-COMPONENT CHAOTIC INFLATION 



In this section we consider a potential with the following form [Lll |2, 



V 



1 " 



2 J.2 

i=l 



(23) 



This potential represents multi-component chaotic infiation. Under slow-roll conditions the masses of the fields are 
similar |2.ll |. and then the number of e-foldings N is given by N{(l),,p) = \ Er=i(<?^i/^p) l2l- The derivatives of N 
are 



r2 



(24) 

N,,=6,,/2M-^. (25) 

Then the contribution in Eq. ^ to /nl is much less than one T^l; {6/5)fNL — ^/N + 0{1/N'^), and therefore it 
is not observable. The higher order contributions in Eqs. H16(l and (|18|l give 



6^/-(i+2) 



16/ 



11 



(n + 2) + : 



(26) 



where r is the tensor to scalar ratio r = 8{H/2Tr)'^ /MpV^, which in this case is given by r = 8/iV. This contribution 
is much smaller than the leading order one in Eq. (jJl, and therefore negligible. 



V. TWO-COMPONENT MODEL 

Here we consider a two-component model given by the following potential 0, ITflj l : 



, 1 </.2(Ar) 1 (j2(7V) 
^+2^^^l^+2^^^lf^ 



(27) 



where r]^ and r]u are the ry slow-roll parameters. 4>{N) and cr(iV) are given in terms of the field values just after 
horizon exit, (p and cr, by 4>{N)_=_6 exp(— Nry^) and o-(iV) = exp(— Nryo-)- Here we consider the case cr = of Ref. 
[Toll. The derivatives of N are IM 



= (lAy^*^), A^a^O, 
= -V4>i^/V4><l>)'^ , 



f?a(l/?700) 



2 2N(»,^->,„) 



(28) 
(29) 
(30) 



The Jnl function is dominated by the term in Eq. (|3J proportional to A^o-o-, and it is approximately given by 



6 



/NL 



16 



(31) 



where e is the slow roll parameter across the a — trajectory, e — {ri^(j)/^/2AIp)^ , and r is the tensor to scalar ratio 
defined earlier. Similarly, A/^^^'' is dominated by the term proportional to iVo-o-, and therefore the main contribution 
comes from A/^]^: 



-Af(^) 



11 / r \3 . 



(32) 



This contribution is much smaller than /nl] 



AfW 
/nl 



-2Ar(r,^-,,„) 



(33) 



VI. CURVATON SCENARIO 



In the curvaton scenario, the curvature perturbation ^ is generated during the oscillations of the curvaton field a 
after the end of inflation |2y, |23, llg ■ The curvaton starts oscillating around the minimum of a quadratic potential, 
V{t,x) = m^(T^(t, x)/2, in an initially radiation dominated epoch with amplitude g = g{a^) until it eventually decays. 
Here tr, is the value of the curvaton field just after horizon exit during infiation. Defining N as the number of e-foldings 
from the beginning of the oscillations to the curvaton decay, its derivatives are given by 

N.,., = ^ [2f7. (2 + n„) g'^ 3 {g'^ + gg")] , (35) 

where g' and g" denote the first and second derivatives with respect to cr* , and is the final fraction of the curvaton 
energy density before it decays. Then 0, Q 

5 5^ 5^ / gg" 
3-6^^+4^- V^-g 



fNL = -^-^^^a + -^A^ + -Za] ■ (36) 



Defining = Ai^ (Vct/\/2F)^, where Va ~ dV/da, the relative contribution of Afj^^^^ to Jnl is 

r P^'/nl , (37) 



'^.(1+2) N 



\ Jnl 

being too small to be detected because of the current observational constraint, Jnl < 121 [T^ . 



VII. CONCLUSION 



We have calculated a higher-order contribution of the three-point correlator of the field perturbations to the three- 
point correlator of the curvature perturbation, which is the next to leading order in the field perturbations. Even 
though intuition suggests that it should be smaller than the leading-order contribution, we have not been able to give 
a general proof that this is so. Instead, we have shown that it is negligible for the cases that have been considered 
so far in the literature. Keeping only the first two terms of Eq.(|21), there is one higher order contribution, and more 
contributions will be generated by higher order terms in Eq.Q. All of these contributions can be calculated in a 
similar way to the calculation shown here. 



Acknowledgments 



IZ is partially supported by Lancaster University Physics Department. YR is supported by COLCIEN- 
CIAS (COLOMBIA). DHL is supported by PPARC grants PPA/G/S/2002/00098, PPA/G/O/2002/00469, 
PPA/Y/S/2002/00272, PPA/V/S/2003/00104 and EU grant MRTN-CT-2004-503369. 



N. Bartolo, S. Matarrese, and A. Riotto, Phys. Rev. D 65, 103505 (2002) arXiv:hep-ph/0112261 . 

V. Acquaviva, N. Bartolo, S. Matarrese, and A. Riotto, Nucl. Phys. B 667, 119 (2003) arXiv:astro-ph/0209156| . 

J. Maldacena, JHEP 0305, 013 (2003) arXiv:astro-ph/0210603l; 

D. H. Lyth and Y. Rodriguez, Phys. Rev. D 71, 123508 (2005) arXiv:astro-ph/0502578''. 
D. H. Lyth and Y. Rodriguez, Phys. Rev. Lett. 95, 121302 (2005) arXiv:astro-ph/0504045 . 
L. Boubekeur and D. H. Lyth, Phys. Rev. D 73, 021301(R) (2006) arXiv:astro-ph/0504046 . 
Y. Rodriguez, PhD Thesis, Lancaster University, Lancaster UK (2005) arXiv:astro-ph/0507701 . 

D. Seery and J. E. Lidsey JCAP 0509, Oil (2005) arXiv:astro-ph/0506056 . 

G. I. Rigopoulos, E. P. S. Shellard, and B. J. W. van Tent, Phys. Rev. D 72, 083507 (2005 ) |arXiv:astro-ph/0410486| . 
G. I. Rigopoulos, E. P. S. Shellard, and B. J. W. van Tent, arXiv:astro-ph/0506704 , arXi v:astro-ph/051104ir 

D. H. Lyth and L Zaballa, JCAP 0510, 005 (2005) arXiv:astro-ph/0507608 . 

L. E. Allen, S. Gupta, and D. Wands, ,]CAP 0601, 006 (2006) arXiv:astro-ph/0509719| . 
L. Alabidi and D. H. Lyth, arXiv:astro-ph/0510441 . 
N. Barnaby and J. M. Cline, arXiv:astro-ph/060148ll 

F. Vernizzi and D. Wands, arXiv:astro-ph/0603799 . 

E. Komatsu et ai, Astrophys. J. Suppl. Ser. 148, 119 (2003) 'arXiv:astro-ph/0302223l. 
M. Sasaki and E. D. Stewart, Prog. Theor. Phys. 95, 71 (1996) arXiv:astro-ph/950700'T]. 
D. H. Lyth, K. A. Malik, and M. Sasaki, .JCAP 0505, 004 (2005) arXiv:astro-ph/0411220]. 

K. Enqvist and A. Vaihkonen, JCAP 0409, 006 (2004) arXiv:hep-ph/0405103 . _^ 

D. Binosi and L. Theussl, Comput. Phys. Commun. 161, 76 (2004) arXiv:hep-ph/0309015 . 
M. Crocce and R. Scoccimarro, Phys. Rev. D 73, 063519 (2006) arXiv:astro-ph/0 509418, . 
S. Dimopoulos, S. Kachru, J. McGreevy and J. G. Wacker, arXiv:hep-th/0507205J. 

D. Polarski and A. A. Starobinsky, Nucl. Phys. B 385, 623 (1992). 
D. H. Lyth and A. Riotto, Phys. Rep. 314, 1 (1999) arXiv:hep-ph/9807278 . 
K. A. Malik, JCAP 0511, 005 (2005) arXiv:astro-ph/0506532 . ' 
D. H. Lyth and D. Wands, Phys. Lett. B 524, 5 (2002) arXiv;hep-ph /Ol 100021 . 

T. Moroi and T. Takahashi, Phys. Lett. B 522, 215 (2001) [Erratum-ibid. B 539, 303 (2002)] arXiv: hep-ph /Ol 10096| . 
D. H. Lyth, C. Ungarelh, and D. Wands, Phys. Rev. D 67, 023503 (2003) arXiv:astro-ph/0208055 . 



